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Abstract. In this paper we investigate a prob- 
lem of state estimation for the dynamical system 
described by the linear operator equation with un- 
known parameters in Hilbert space. We present 
explicit expressions for linear minimax estima- 
tion and error provided that any pair of uncer- 
tain parameters belongs to the quadratic bound- 
ing set. As an application of the main result we 
present the solution of minimax estimation prob- 
lem for the linear descriptor differential equation 
with constant matrices. 
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Introduction 

One of the major problems in applied mathe- 
matics is the problem of state estimation for dy- 
namical systems described by the linear equations 
with uncertain parameters. This problem belongs 
to the so called "uncertain inverse problems". In 
the strict sense the problem is described as fol- 
lows: given some element (for instance measure- 
ments of the system output) y from some func- 
tional space one needs to estimate the expression 
£(9) provided that 9 obey the equation g{9) = 0. 
This problem is non-trivial if there exists non- 
unique 9 satisfying the equation g(9) = and 
the equality y = C{9) holds. In this case the 
estimation problem may be reformulated as fol- 
lows: given y = C(9),9 G Q,y G Y one needs 
to find the estimation £(9) of the expression £{9) 
provided that g(9) = and C(-),£(-) are given 
functions. Note that this problem becomes trivial 



if the equation y = C(9) has the unique solution 

9. Really, in this case we just set £{9) := £{9). 

The estimation problem is said to be linear if 
0,1" are linear spaces and C(-),£(-) are linear 
mappings. It is the common case when 

C{9) =H<p + Dr,, g{6) = Lp + Bf, (*) 

where 9 = (x, f, rj) C X, F, Y, H, D, L, B are lin- 
ear mappings. The linear estimation problem is 
said to be uncertain if D ^ 0, L and B are non- 
trivial or if B = then N(L) = {ip : Lip = 0} ^ 
{0}. Note that the choice of solution method de- 
pends on the "type of uncertainty" : if /, n denotes 
realizations of random elements then it's natural 
to apply probability methods. This requires an 
a priori knowledge of distribution characteristics 
of the random elements. In the sequel we assume 
that there is uncertainty in (*) if distributions of 
random elements or some deterministic parame- 
ters of the system are partially unknown. An up 
to date description of the state of the art in the 
theory of uncertain estimation problems with spe- 
cial £, L, H, B, D in special spaces is to be found 

at eh ei sua ei e]. 

Classical theory of uncertain estimation prob- 
lems [H El S] works well when the linear mapping 
L in (*) has bounded inverse mapping. The solu- 
tion of the linear uncertain estimation problem for 
linear equations with Noether operator was intro- 
duced in [7j. Note that the introduced approach is 
not suitable when dim(L) = oo. This is the case 
for linear differential-algebraic (DAE) equations 

0- 

The major contribution of this paper is an state 
estimation approach for uncertain equations with 
linear closed operator in Hilbert space. This ap- 
proach is valid when L has non-closed range. It 
generalizes the classical theory [H El H] to the 
linear noncausal uncertain differential-algebraic 
equations. The application of presented method 
to the linear equations with Noether operator was 
described in [5]. 

Notation. Set c(G, 
G}, let S(G, •) denotes 
of Q, set dom/ = {x 
f*(x*) = sup x {(x*,x) 
inf{c(G, z), L*z = u}, 



= sup{(z,/),/ G 
the indicator function 
G H : f(x) < oo}, 
- /(*)}, (L*c)(u) = 
(fL)(x) = f(Lx), 



1 to appear in Ukrainian mathematical journal, 2009 (L*c)(u) = inf {c(G, z) , L* Z = u], let c\f = f* 
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denotes the closure of / in the Young-Fenhel 
sense, Arginf u /(w) denotes the set of minimum 
points of /, Pi* denotes the orthogonal projector 
onto R(L*), df(x) denotes the subdifferential of / 
at x and (•, •) denotes the inner product in Hilbert 
space. 

Problem formulation and def- 
initions. 

Suppose that Lp G £f and 

y = Hip + rj (1) 

The mappings L, H and the set £f are supposed 
to be given. The element 77 is uncertain. Our 
aim is to solve the inverse problem: to construct 
the operator mapping the given y into the estima- 
tion l{(f) of expression £((p) and to calculate the 
estimation error a. Now let us introduce some 
definitions. 

The operator L : H i— > T is assumed to be 
closed. Its domain 3>{V) is supposed to be a dense 
subset of the Hilbert space TC, H G Jf(Tl,y). 
Note that the condition Lp G £f is equal to the 
following 

L<p = f, (2) 

where / is uncertain and belongs to the given sub- 
set £f of the Hilbert space T . In the sequel 77 is 
supposed to be a random 3^-valued vector with 
zero mean so that its correlation R v G ffl, where 
M is some subset ofJf(y,y). Also we deal with 
deterministic 77 so that (/, 77) G Q, where Q is some 
subset of T x y. Note that the realization of 
y depends 77, H and /. Also it depends on ele- 
ments of N(L) = {ip G 9{L) : Lp = 0} so that 
y = H((po + <p) + r), where <po may be thought as 
inner noise in the state model (J2J). 

Let £((p) = (£, (p), £(<p) = (u, y) + c. Since L, H 
are not supposed to have a bounded inverse map- 
pings the £(<p>) and £{<p) are not stable with a re- 
spect to small deviations in /, rj. Also /, rj are 
supposed to be uncertain. Therefore we use the 
minimax design in order to construct the estima- 
tion. 

Definition 1. The function £((p) = (u, •) + c is 
called the a priori minimax mean-squared estima- 



tion iff a(£,u) = inf UiC a(£, u) where 

a(£,u):= sup M(%) -%)) 2 (3) 

The number cr(£) = a^(£, u) is said to be the min- 
imax mean-squared error in the direction £. 

On the contrast the a posteriori estimation de- 
scribes the evolution of the central point of the 
system reachability set 

(Lcp,y- Htp) G Q 

consistent with measured output y [TQl [TTJ, [T2] . 
Note that the condition (Lip, y — Hip) G Q holds 
if Hy 1 1 < C for some real C. But it doesn't hold 
in our assumptions if 77 is random since ||-R^|| < 
c doesn't imply ||y|| < C for realizations of rj. 
Therefore rj is supposed to be deterministic. 

Definition 2. The set 

X y = {p G @{L) : {Lip, y-Hp)E G} 

is called an a posteriori set. The vector ip is said 
to be minimax a posteriori estimation of ip in the 
direction I iff 

d(t) := inf sup \(£,p)-(£,ip)\ = sup \(l,<p)-(£,i>)\ 

The expression d(£) is called the minimax a pos- 
teriori error in the direction £. 

In the sequel the minimax mean-squared a pri- 
ori estimation (error) is referred as minimax esti- 
mation (error). 

Main results 

Proposition 1. Assume that M are convex 
bounded closed subsets of T , JZ'iy, y) respec- 
tively. For the given £ G Ti the minimax error 
a(£) is finite iff 

£ - H*u G dom c\(L*c) fl (-l)dom d(L*c) (4) 

for some u G y . Under this condition 

cr(£,u) = sup (R v u,u)-\- 

-[cl(L*c) (£ - H*u) + c\(L*c) (-£ + H*u)} 2 
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where Corollary 2. Assume that linear mappings L : 

U^T, H elf(H,y) obeyl) or 2) (Cor.ty. 
\ ) c domcl(L c) C K{L ) Then p has fl soton | e 9(L*),p G ^(L) 

£ = L*z + /or some z G 3>(L*),u G 3^ 
Corollary 3. Under the conditions of Cor. U\for 



If Aigmf u a (£,u) ^ 0, then £((p) = (u,y) + c, 
where 



and 



a » . r //, \ any £ G R(L*) + R(H*) and some realization of 

u>e aave = \£,(p), where (p obeys the 



equation 



c=^{cl{L*c){£-H*u)-cl{L*c){-£ + H*u)) L*q = H*(y - H0), 

L0 = q 



Consider an a posteriori estimation. 



Theorema 1. Suppose that&i is convex bounded 

closed balanced set and its G int Also assume 

that Proposition 2. Let Q be a convex closed bounded 

7] G {r] : M(r],r]) < 1} subset ofyxj 7 . The minimax a posteriori error 

„, , „ . „ ^ , ., . . .. in the direction £ is finite iff £ G domc(A'„, •) fl 

1 hen for the given t G n the minimax estima- , . N . , . . N , v " ' 

«//)\ • n -re n tt* n/r*\ r ( — 1 IQOHI C[A. V , • ) and 

tion a{£) is finite iff £ - H*u G R(L*) for some v ' v v ' 

u G y . Under this condition there exists a unique R(L*) + R(H*) C domc(X y , •) fl (— l)domc(X y , ■) C 

minimax estimation u <EUg and R(L*) + R(H*) 

a(£,u) = mma(e,u), (10) 

u 

a(£, u) = (u, u) + min{c 2 (^i, z),L*z = £- H*u} Under this condition 

(6) (£,<£) = ~(c(X y ,e)-c(X y ,-£)), 

If R{L) , H (N (L)) are closed sets then u is deter- _ z( c (X £) + c(X —£)) 

mined by the following conditions 2 

Theorema 2. Let 

u-Hp e H(dI 2 (H*u)),Lp = 0, 

I 2 (w) = mm{c 2 (& u z),L*z = P L *(£-w)}, (7) G = {(f,v) ■ W + IMI* < 1}, 

z 

and assume that one of the conditions of Cor. [I] 

Corollary 1. Let ^ s f u ifiH e d, 77^ minimax a posteriori estimation 



^ = {/ G T: (/,/) < 1},77 G {7/ : M(^) < 1}, 



<£> o&eys i/ie equation 



J thf L*q = H*(y-H0), 

and suppose that r „ 

Ly? = 

i; R(L),H(N(L)) are closed sets; . ff £ g ^ + ^ ^ error , s 

i?(T) = {[La;, ifx] , x G @(L)} is closed set. 9 wen h V 

Then only for £ G R(L*) + R(H*) the unique mm- d{l) = (1 - (y,y - H<p))ia(£) (13) 

imax estimation u is given by u = Hp, where p is Corollary 4. Assume that the conditions of The- 

any solution of the equations orem m are fulfilM and ^ = ^ ^ for my ^ 

L*2 _ £ _ u*}jp where (p obeys ffT^l) . Then <p gives the minimax a 

T ^ ^ (8) posteriori estimation of so that 

Lp = z ^ 

inf sup \\<f — x\\ = 

The minimax error is given by the following ex- ^ x v xex y 

pression sup ||(^ — x|| = (1 — (y, y — if(^)^ max a(£) 

a(£) = (£,p)z x€X y IKII=i 
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Now we shall apply Cor. H]to the linear uncer- 
tain estimation problem for differential-algebraic 
equation. Taking into account that any DAE with 
constant matrices has the SVD canonical form [13] 
we assume without loss of generality that 

r — \ o oj> b - lc 3 C A ) 

Proposition 3. Lett \— > x(t) € M. n be the solution 
of 



Fx(t)-Cx(t) = f(t),Fx{to) = 0, 



d_ 

It 

and set 



e = {{f,v)-- / (ll/(*)lr + N(*)lr)*<i} 

J to 

Then the minimax a posteriori estimation of the 
function x(-) is given by x(-), where x(t) = 
[xi(t),x 2 {t)\, 

Xl (t) = (Ci - c 2 (e + c 4 c 4 y 1 c 4 c 3 )x 1 (t)+ 

(C 2 (E + C' 4 C 4 y 1 C' 2 + E)q 1 (t) + 
C 2 (E + C' 4 C,)- l y 2 (t), 

qi (t) = {-C[ + C' 3 C 4 (E + C' 4 C 4 )- 1 C' 2 )q 1 (t) + 

C'. i C 4 {E + C' 4 C 4 y 1 y 2 {t)-y l {t) + 

(C' 3 (E - C A {E + C^r^Ca + E)x 1 (t), 

q 1 (T) = 0,x 1 (t ) = 0, 

x 2 (t) = —{E + C' 4 C 4 )- 1 C' 4 C 3 x 1 (t)+ 

(E + C;C 4 )- 1 (C^ 1 (t)+y 2 (t)), 

q 2 (t) = —(E - C 4 (E + C;C 4 )- 1 C;)C 3 x 1 (t)- 



C 4 {E + C' 4 C 4 )-\C' 2qi {t) + y 2 {t)) 



(14) 



provided that y(t) = x(t) + rj(t), t < t < T . The 
minimax error is given by 

sup \\x— x\\ = (1— / (y, y— x)di) a max( / (£,p)dt) 
x y Jto n =l Jt 

where p(-) obeys f|T4]) with y(t) = £(t). 
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